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Abstract 

The restricted circular three-body problem is considered for the fol- 
lowing parameter values C — 3.03, = 0.0009537 - the values for Oterma 
comet in the Sun- Jupiter system. We present a computer assisted proof 
of an existence of homo- and heteroclinic cycle between two Lyapunov 
orbits and an existence of symbolic dynamics on four symbols built on 
this cycle. 



1 Introduction and statement of results 

In paper jll| methods of dynamical system theory were used (see also |l3|) to 
explain rapid transitions from heliocentric orbits outside the orbit of Jupiter to 
heliocentric orbits inside the orbit of Jupiter and vice versa for Jupiter comets 
Oterma and Gehrels 3. To model this problem authors in JTTI used planar 
circular restricted three-body problem and established that for a parameters 
corresponding to Sun-Jupiter-Oterma system rapid transitions of Oterma are 
explained by transversal intersections of stable and unstable manifolds of two 
periodic orbits around libration points Li and L2- In fact the existence of 
symbolic dynamics on three symbols was claimed. 

The goal of this paper is develop and test tools which allow with computer 
assistance to prove the results claimed in [ pT| . 

Before we state our main results we give a short description of the planar 
restricted circular three-body problem. We follow the paper and use the 
notation introduced there. Let S and J be two bodies called Sun and Jupiter, 
of masses TOs = 1 — and rrij = fi, ^ G (0, 1), respectively. They rotate in the 
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plane in circles counter clockwise about their common center and with angular 
velocity normalized as one. Choose a rotating coordinate system (synodical 
coordinates) so that origin is at the center of mass and the Sun and Jupiter 
are fixed on the x-axis at {—^J., 0) and (1 — /z, 0) respectively. In this coordinate 
frame the equations of motion of a massless particle called the comet or the 
spacecraft under the gravitational action of Sun and Jupiter are (see and 
references given there) 



where 



X -2y = nx{x,y), y + 2x ^ ny{x,y), (1.1) 



n{x,y) 



n r2 



n = \J {x + /i)2 + r2 = v/(2: - 1 + m)^ + 



y 



Equations (1.1) are called the equations of the planar circular restricted three- 
body problem (PCR3BP). They have a first integral called the Jacobi integral, 
which is given by 

C(x, y, i, y) = -(±2 + y2) + 2r!(x, y). (1.2) 

We consider PCR3BP on the hypersurface 

M{^i, C) = {(x, y, X, y) \ C(x, y, x, y) = C}, (1.3) 

and we restrict our attention to the following parameter values C = 3.03, fi = 
0.0009537 - the parameter values for Oterma comet in the Sun- Jupiter system 



(see fn)). 

The projection of A^(/i, C) onto position space is called a Hill's region and 
gives the region in the (x, ?/)-planc, where the comet is free to move. The Hill's 
region for the parameter considered in this paper is shown on Figure |l| in white, 
the forbidden region is shaded. The Hill's region consists of three regions: an 
interior (Sun) region, an exterior region and Jupiter region. 

In |ll| a very good numerical evidence was given for the following facts for 
the Sun-Jupiter-Oterma system: 

0. an existence of Lyapunov orbits and L2 around libration points Li and 

L2, respectively. Both orbits arc hyperbolic and are located in the Jupiter 
region. 

1. There exists a transversal heteroclinic orbit connecting LI and L^. There 

exists a transversal heteroclinic orbit connecting L2 and LJ. Both orbits 
are in the Jupiter region. These orbits were discovered for the first time 
in0. 

2. there exists a transversal homoclinic orbit to in the interior (Sun) region, 

3. there exists a transversal homoclinic orbit to in the exterior region. 
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Figure 1: Hills region for PCR3BP with C = 3.03, /i = 0.0009537 frora |ri|. 



By transversal hetero- and homoclinic orbit, we mean that appropriate unstable 
and stable manifolds intersect transversally. For example in assertion 1: the 
stable manifold of L\ intersect transversally the unstable manifold of ■ 

It is now standard in dynamical system theory (see |0 and references given 
there) to derive from assertions [0-3] an existence of symbolic dynamics on four 
symbols S,L\,L2,X which the following allowed transition 

In pl|( section 1.4 ) an existence of symbolic dynamics on three symbols, only, 
was claimed. Instead of two symbols L\ and L\ one symbol J for Jupiter region 
was used. 

From the point of view of rapid transition of Oterma from interior region 
to an exterior region and vice versa, an existence of heteroclinic orbits between 
L\ and L\ claimed in assertion 1 was of a special importance, as they are an 
indication of an existence of a dynamical channel joining an interior with an 
exterior region. 

The following two theorems summarize main results of our paper 

Theorem 1.1 For PCR3BP with C = 3.03, = 0.0009537 there exist two 
periodic solutions in the Jupiter region, L\ and L2 , called Lyapunov orbits, and 
there exists heteroclinic connections between them, in both directions. Moreover 
for both orbits L\ and L2 there exists a homoclinic orbit in interior and exterior 
region, respectively. 

The next theorem says that the homo- and heteroclinic connections whose 



existence is established in Theorem 1.1 are topologically transversal i.e. give 



rise to the symbolic dynamics, just as in the case of an existence of transversal 
intersections of stable and unstable manifolds. 

Theorem 1.2 For PCR3BP with C 3.03, /i = 0.0009537 there exist a sym- 
bolic dynamics on four symbols {S, X, Li, L2} corresponding to Sun and exterior 
regions and vicinity of Li and L2, respectively. 
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A precise statement of this theorem with aU necessary details about the symbohc 
dynamics is given as Theorem 7.1 in Section 0. 

Hence we have proved assertion 0, but we did'nt proved assertions 1-3, as 
we did'nt checked that stable and stable manifolds of Lyapunov orbits intersect 
transversally. Instead we had proved that there is enough topological transver- 
sality present to build a symbolic dynamics on it. The use of topological tools 
was essential for the success of this work, as the rigorous computation of stable 
and unstable manifolds appears to be much more difficult than the computations 
reported in this paper. 

Contents of our paper may be described as follows. In Section ^ we continue 
our brief description of PCR3BP, we define suitable Poincare maps and state 
their symmetry properties. In Section || we present the main topological tool 
used in this paper - the notion of covering relation. In Section Q we describe 
how to link a local hyperbolic behavior with covering relations to obtain homo- 
and heteroclinic orbits. In Sections ^ and || we repor t th e results of our rigorous 
computations for PCR3BP and we prove Theorem |l.l[ In Section ^ we show 
how to use symmetries of PCR3BP together with covering relations to complete 
the proof of Theorem 1.2. Section ^ contains the details of the numerical part of 
proof, we mainly discuss the question of an efficient approach to a verification of 
covering relations. We also include all initial data, so that a willing reader with 
his own code can verify our claims. In Section ^ we provide a brief description 
of basic classes used in our programm and we give the names of functions, which 
perform proofs of numerical lemmas, in order to facilitate an understanding of 
the source code, which we had made available on-fine (see||l^). In Section ^ 
we discuss some natural extensions of our results. 

What is new in this paper besides giving a proof for some results from Jill? 
First of all it shows how to successfully link numerically cheap C°-methods (the 
covering relations) with much more numerically expensive C-'^-methods (a local 
hyperbolicity). This was previously done for maps (see |^), only. The main 
obstacle in applications to ODEs was the lack of an efficient ODE solver. 
Such a solver - called a C^-Lohner algorithm - was recently proposed by the 
second author in |2|] . 

Other novelties in this paper are some theoretical improvements in the theory 
of covering relations. We use an abstract definition of covering relation from 
[0 and we show how to use a symmetry which involves a time reflection (for 
Poincare map this correspond to taking an inverse map) together with covering 
relations. Both these improvements (in numerical algorithms and in theory of 
covering relations) allow to reduce the computation time considerably (probably 
by two orders of magnitude). The total computation time on a PC using a 
l.lGhz Celeron processor is less than 40 minutes. 
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2 Properties of PCR3BP: Poincare maps and 
symmetries 

In this section we continue our brief description of PCR3BP we started in the 
Introduction and we introduce various notations which wiU be used throughout 
the paper. 

The PCR3BP has three unstable colhnear equihbrium points on the Sun- 
Jupiter hne, caUed Li, L2 and L3 (see Fig. 1.4 in ), whose eigenvalues 
include one real and one imaginary pair. The value of C (Jacobi constant) at 
the point Li will be denoted by Ci. An linearization at Li for i — 1^2 for 
the parameter range considered here, shows that these points are of center- 
saddle type (see jl^). By theorem of Moser jl^ it follows that for C < Ci and 
|C — Ci| small enough, there exist hyperbolic periodic orbits, L*, around Li, 
called Lyapunov orbits. Observe that for a fixed value of C < Ci an existence 
of Lyapunov orbit L* is not settled by the Moser theorem and has to be proved. 

As was mentioned in the Introduction we restrict our attention to the fol- 
lowing parameter values C — 3.03, /i = 0.0009537 - the parameter values for 
Oterma comet in the Sun- Jupiter system (see |]TT|). Since we work with fixed 
parameter values we usually drop the dependence of various objects defined 
throughout the paper on /x and C, so for example Ai = A4{fi,C). For our 
parameter values we have C2 > C > C3 ( this means that we are considering 
Case 3 from section 3.1 in ||TT[|). 

We consider Poincare sections: Q — {{x,y,x,y) G M\y = 0}, 9+ = On{?) > 

0}, e_ = en{y<o}. 

On Q± we can express y in terms of x and x as follows 

y^±y/2n{x,0)-x^ (2.4) 

Hence the sections 0± can be parameterized by two coordinates {x, x) and we 
will use this identification throughout the paper. More formally, we have the 
transformation T± : ^ 0± given by the following formula 

T±{x,x) = {x,0,x,±y/2n{x,0) - x^ - C ) (2.5) 

The domain of T± is given by an inequality 2n{x, 0) — x^ — C > 0. 

Let Tr± : Q± — > R and tt^ : Q± — > M denote the projection onto x and x 
coordinate, respectively. We have Tr±{xQ,XQ) = xg and iTxixo^xo) = xq. 

We will say that {x, x) G 0± meaning that {x, x) represents two-dimensional 
coordinates of a point on 0±. Analogously we give a meaning to the statement 
M C e± for a set M C M^. 

We define the following Poincare maps between sections 

p+ : e+ ^ e+ 
p_ : e_ ^ e_ 
Pi,+ : e+ ^ 0_ 

Pi _ : 0_ ^ 0+. 
2 ' 
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As a rule the sign + or — tells that the domain of the maps P± or Pi _|_ is 
contained in Q± (the same sign). Observe that 

P+{x) = o P.^^ix), P^{x) = Pi,+ o P. _{x) 

whenever P+{x) and P-{x) are defined. These identities express the following 
simple fact: to return to 0+ we need to cross Q with negative y (this is Pi _|_ 
first and then we return to O with y > (this is Pi _). 

Sometimes we will drop signs in P± and P^ _|_, hence P{z) — P+{z) if z G 8+ 
and P{z) = P-{z) if z e a similar convention will be applied to Pi. 

2.1 Symmetry properties of PCR3BP 

Notice that PCR3BP has the following symmetry 

R{x, y, X, y, t) = {x, -y, -i, y, -t), (2.6) 

which expresses the following fact, if {x{t),y{t)) is a trajectory for PCR3BP, 
then (x(— t), —y{—t)) is also a trajectory for PCR3BP. From this it follows 
immediately that 

if P±{xq,Xo) = {xi,ii), then P±{xi, -xi) = {xq, -xq) (2.7) 
if Pi ±(a;o,io) = (a;i,a;'i), then Pi -^{xi,~xi)^{xo,-xo) 

We will denote also by R the map R : 8± Q± R{x, x) = {x, —x) for 
{x,x) e 6±. Now eq. ( ^!^ ) can be written as 

if P±(a;o) = a;i, then P±{R{xi)) = R{xq) (2.8) 
if Pi^±(a:o) = xi, then Pi -j-(P(xi)) = R{x()) 

3 Topological tools 

In this section we present main topological tools used in this paper. The crucial 
notion is that of a covering relation. This notion in various forms was introduced 
in papers ||2^, |2l|, |2^, |2^ . Here we follow the most recent and most general 
version introduced in [[7| and the reader is referred there for proofs. 

3.1 h-sets 

Notation: For a given norm in R" by P„(c,r) we will denote an open ball of 
radius r centered at c G M" . When the dimension n is obvious from the context 
we will drop the subscript n. Let 5'"(c,r) = 9P„+i(c, r), by the symbol S*" we 
win denote S'"(0, 1). We set R[_= {0}, Po(0, r) = {0}, dBo{0, r) = 0. 

For a given set Z, by int Z, Z, dZ we denote the interior, the closure and the 
boundary of Z, respectively. For the map h : [0, 1] x Z — + M" we set ht = h{t, •). 
By Id we denote an identity map. For a map /, by dom(/) we will denote 
the domain of /. Let / : C M" M" a continuous map we will say that 
X C dom (/^^) if the map /^^ : X M" is well defined and continuous. 
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Definition 3.1 A h-set, N, is the object consisting of the following data 

• \N\ - a compact subset ofW^, a support of N 

• u{N), s{N) £ {0, 1,2,...}, such that u{N) + s{N) = n 

• a homeomorphism cn ■ K" ^ M" = ]R"(^) x R<'^\ such that 

CNi\N\) = :b;:^(0, 1) X B:^(0, 1). 

We set 

iVc = B;:^(o,i) x;b:^(o,i), 

= aB;:^(o, 1) X 51^(0, 1) 
N+ = 'K^){o, 1) X as:^(o, 1) 

N-=c],\N-), N+=c],'{N+) 

Hence a h-set, N, is a product of two closed balls in some coordinate system. 
The numbers, u(N) and s{N), stand for the dimensions of nominally unstable 
and stable directions, respectively. The subscript c refers to the new coordinates 
given by homeomorphism cjy. We will call (N^) an exit set of N and 
{N+) an entry set of N. Observe that if u{N) = 0, then A^- = and if 
s{N) = 0, then N+ = 0. 

Definition 3.2 Let N be a h-set. We define a h-set N'^ as follows 

• |iV^| = |iV| 

• u{N'^) = s{N), s{N^) = u{N) 

• We define a homeomorphism c^t : M" -> R" = M"(^^) x M''(^'^), by 

Cmt{x) ^j{cN{x)), 

where j : R"(^) x W'-^'^ R"(^) x R"(^) ts given by j{p, q) = {q,p). 

I 

Observe that A^-^'+ = N~ and N'^'~ = N~^. This operation is useful in the 
context of inverse maps, as it was first pointed out in 

3.2 Covering relations 

For n > and a continuous map / : S'" ^ 5" by d{f) we denote the degree 
of / [Q. For n = we define the degree, d{f), as follows. Observe first that 
S° = {-1,1}. We set 

fl, if /(I) = land /(-I) = -1, 
dif)^l-l, if/(l)==-land/(-l) = l, (3.9) 
[O, otherwise. 
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Definition 3.3 Assume n > 0. Let f : B„(0, 1) M", such that 

0^/(95(0,1))). (3.10) 



We define a map Sf : ^ — > S"" ^ by 



Definition 3.4 Assume N, M are h-sets, such that u{N) = u{M) = u and 
s{N) = s{M) — s. Let f : \N\ R" be continuous. Let = cm ° f ° c'^ : 
A^c R" X R" . Let w be a nonzero integer. We say that 

(N f -covers M with degree w) iff the following conditions are satisfied 

1. there exists a continuous homotopy h : [0, 1] x Nc —> R" x R^, such that the 
following conditions hold 

ho = fc, (3.12) 
h{[0,l],N-)nMc = (3.13) 
h{[0,l],Nc)nM+ = (3.14) 

2.1 If u> 0, then there exists a map A : M" R", such that 

hi{p,q) = {A{p),0), where p^W and q(^W (3.15) 
A(9B„(0,1)) C R" \ S„(0, 1) (3.16) 

Moreover, we require that 

d{sA) = 

2.2 Ifu = 0, then 

hi{x) = forxeN^ (3.17) 
w = 1. (3.18) 

Intuitively, N => M if / stretches N in the 'nominally unstable' direction, 
so that its projection onto 'unstable' direction in M covers in topologically 
nontrivial manner projection of M. In the 'nominally stable' direction N is 
contracted by /. As a result N is mapped across M in the unstable direction, 
without touching M"*". An example of covering relation on the plane with one 
unstable direction is shown on Figure |[ 
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Definition 3.5 Assume N, M are h-sets, such that u{N) = u{M) — u and 
s{N) = s(M) = s. Let g -.new ^W^. Assume that g-^ : \M\ R" is well 
defined and continuous. We say that N ^= M (N g-backcovers M with degree 

w) iffMT& N^. 

The following theorem proved in is one of main tools used in this paper. 
Various versions of this theorem (without backcovering) using slightly weaker 
notions of covering relations or even without an explicitly defined notion of 
covering relation were given in 1^, ^ p^ , ||. In the planar case this theorem 
with backcovering was stated also in [ ij. 

Theorem 3.6 Assume Ni, i = 0, . . . , fc, Nk — Nq are h-sets and for each 
i = 1, . . . ,k we have either 

iV,_i ^ (3.19) 



or \Ni\ C dom{f^ ^) and 



iV,„i N,. (3.20) 



Then there exists a point x G int |iVo|, such that 

/zo/,_io...o/i(a;) e mt\N,l i = l,...,k (3.21) 
fkofk-io---ofi{x) = X (3.22) 

Obviously we cannot make any claim about the uniqueness of x in Theo- 
rem B.6, 



3.3 Covering relation on the plane with one nominally 
expanding direction (tt = 1) 

In this section we discuss the case, when u — s — 1, hence we have only one 
nominally expanding and one nominally contracting direction. The basic idea 
here is: the set N~ consists from two disjoint components and all possible 
values of the degree w in covering relation are ±1. This allows to give sufficient 
conditions for an existence of covering relations, which are relatively easy to 
verify. 

Definition 3.7 Let N be a h-set, such that u{N) = s{N) — 1. We set 

Ni- - {-l}x[-l,l] 

- {l}x[-l,l] 

5(iV)^ = (-CX), -l)xR 

S{N)l = (l,cx))xM. 



9 



Figure 2: An example of h-set on the plane. 



We define 



We will call iV'^ S{Ny and S{NY the left edge, the right edge, the left 

side and right side of N , respectively. 

It is easy to see that N' = N ^" U iV^ 

The triple {\N\,S{Ny,S{Ny) is a t-set from [|. As in |] we will use the 
following notation for S{NY''-. 



Remark 3.8 For all h-sets used in this paper the support is a parallelogram. A 
usual picture of a h-set is given in Figure |^. 

A typical picture illustrating covering relation on the plane with one 'unsta- 
ble ' direction is given on Figure ^. 

The following theorem was proved in |0] for any n > 1 and u{N) = 1. Here 
we rewrite it for the planar case in a slightly different notation (we use A^' and 
N"- for S{Ny and S{MY, respectively). 

Theorem 3.9 Let n = 2 and let N, M be two h-sets in M", such that u{N) = 
u{M) = 1 and s{N) = s{M) = 1. Let f : \N\ R" be continuous. 

Assume that there exists qq G Bs{0,1), such that following conditions are 
satisfied 



S{NY, 



N'' = S{NY 



/(cAr([-l,l] X {go})) C int(Af'u|M|U7ir) 
/(|A^|)nM+ = 0, 



(3.23) 
(3.24) 
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Figure 3: An example of an /—covering relation: A^o =^ Nq, Ni 



and one of the following two conditions holds 



/(7V'^) C A/' and f{N'"') C W (3.25) 
/(7V''=) C A'f and f{N''^) C M'. (3.26) 



Then there exists w = ±1, such that 

3.4 Representation of the h-sets 

In this paper we use very simple h-sets, namely the support is a parallelogram. 

A h-set is defined by specifying the triple N = t{c, u, s), where c,u,s € R^, 
such that M, s are linearly independent. We set 

|7V| = {x eR^ \3t^,t,e[-l.l] x^c + tiu + t2s} 

= C+ [-1,1] -11+ [-1,1] -s 
CN{tl,t2) = C + tiU + t2S 

In this representation c is a center point of the parallelogram N , u represents an 
oriented half-length in the 'unstable' direction and s is an oriented half-length in 
the 'stable' direction. See Fig. ^ for an example of h-set in this representation. 
We have 

N^'' = c-u+[-l,l]-s 
N''" = c + M+[-l,l]-s 
iV' — c + (— oo, — l]u + (— oo, cx))s 
iV" — c + [1, oo)u + (— cxD, cxd)s. 
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Figure 4: An example of an i?-symnietric h-set 



We introduce notions of top and bottom edges of N, N and N by 

iV'"' = c+[-l,l]-u-s 
N^" = c+[-l,l]-u + s 



Let us recall that the symmetry R : M.^ 
given by 



introduced in section 2A was 
R{xi,X2) = (xi, -2:2) 



Definition 3.10 A h-set, N, will be called an i?-symmetric h-set ifN = t{c, u, s) 
for some c,u, s G R^, such that 

R{c) = c 
R{u) = s or R[u) = — s 

Figure ^ shows an example of a i?-symmetric h-set. Symmetry properties of 
such h-set are apparent. 

3.5 Action of R on h-sets 

The symmetry of Pi/2,± and P± expressed in ( ^.7|) relates the maps and theirs 
inverse, hence beside mapping the support of TV by i? it will switch also the 
nominally stable and unstable directions. This motivates the following definition 
of the action of the symmetry R on h-sets 

Definition 3.11 Let N be a h-set. We define a h-set R{N) as follows 
. \RiN)\^Ri\N\) 

• u{R{N)) = s{N) and s{R{N)) = u{N) 

• the homeomorphism C7?(7v) : R" ^ ^^W^)) x M''('S'(^)) is given by 

cr{n) = cjvT o R^^ 
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Observe that according to above definition we have 

R{N)- ^ R{N+) = R{N^'-) 
R{N)+ = R{N-) = R{N'^-+) 
R{t{c, u, s)) = t{R{c),R{s),R{u)) (3.27) 



As an immediate consequence of equation (3.27) we obtain 

Lemma 3.12 Let N — t{c, u, s) be an R-symmetric h-set. Then R{N) = N. 

We have the foUowing easy lemma 

Lemma 3.13 Let fi : dom(fi) C ^ for i — 1,2 continuous and invert- 
ible on some open sets. Assume that 

if fi{x)=xi, then f2{Rxi) = R{x). (3.28) 

Let Ni, N2 be h-sets such that 



Then 



Ni N2. 



R{N2) 4^ R{Ni) 



Proof: From Dcf. 3.11 and the assumed symmetry it foUows immediately that 



R{Nif^R{N^). (3.29) 



From above Lemma and ( p. 8]) we obtain 

Corollary 3.14 If Ni ^ N2, then R{N2) ^ R{Ni). 

-Pi + Pi _ 

IfNi ^ N2, then R{N2) 4= R{Ni). 

-Pi _ Pi + 

IfNi ^ N2, then R{N2) <^ i?(iVi). 

4 tools 

The goal of this section is to describe the tools which allow in the presence of 
hyperbolic fixed points for a map to prove an existence of homo- and heteroclinic 
trajectories. 

In this section we recall the results from [|j with some additions (see also 
where the method was outlined for the first time) . In the symbol of covering 

relation we will drop the degree part, hence we will use N =^ M instead of 

N M for some nonzero w. 
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4.1 General theorems 

Let P : M" M" be a C^-map. For any set X we define an interval matrix 
[DP{X)] C M"^" to be an interval enclosure of DP{X) given by 

M e [DP{X)] iff inf DP{x),.j < M,j < sup DP{x)^j i, j = 1, 2, . . . , n 



Lemma 4.1 Let N be a convex set. Assume Xq,xi e N. Then 

- P{xo) e [DP{N)] ■ (xi - xo). (4.30) 
Moreover, there exists a matrix M £ [DP{N)] such that 

Pixi) - Pixo) = M ■ (xi - xo) (4.31) 

Proof: 

^1 dP. 



f dP 

P{xi)-P{xo)= —{xo + t{xi-xo))dt 
r-i Qp 



, (xo + i(a;i - a;o))c^i • - a;o) 
c*^ 



To finish the proof observe that 

^ dP, 



J -g^ixo + t{xi - xo))dt e [DPiN)] 

Let /d : M" ^ R" denote the identity map. 
Theorem 4.2 Let N be a convex set. Assume that 

i det{[DP{N)] - Id) = {t = det{M - Id) \ M e [DP{N)]}, 
then N contains at most one fixed point of P. 

Proof: Assume that xo,xi N are fixed points for P. Then from lemma 
follows that 



4.1 



it 



Xi-Xo = Pixi) - P{xo) = M ■ {xi - Xo) (4.32) 

for some matrix M € [DP{N)]. Hence xi — xo is in the kernel of M — Id. From 
our assumption it follows that det(Af — Id) ^ 0, hence xi = xo. ■ 

Consider a two-dimensional function f{x) ~ {fi{x), f2{x))'^ , where x = 
{xi,X2)'^- We assume that /(O) = 0, i.e. is a fixed point of /. For a convex 
set U, such that G ?7 we define intervals Xi{U), ei{U), £2{U) and X2{U) by 
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Since /(O) = then from Lemma [4.1| it follows that 

e Ai(;7)xi+£i({/)a;2 

f2{x) e e2{U)xi + X2{U)X2 

Let 

e[{U) = sup{|£| : e G £i(C/)}, e'^iU) = snp{\e\ : e G £2(t/)}, 
A;(C/) = inf{|Ai| : Ai e Ai(C/)}, A^(C/) = sup{|A2| : A2 G A2(C/)}. 

Let us define the rectangle -/Vq^ q.^ by 

Definition 4.3 Def. 1] Let be a fixed point for the map f . We say that 
f is hyperbolic on N 3 x^,, if there exists a local coordinate system on N , such 
that in this coordinate system 

= (4.34) 
e[{N)e'2{N) < {1 - X',{N)){X[{N) - 1). (4.35) 
N = iV„,,„,, (4.36) 

where ai > 0, q;2 > are such that the following conditions are satisfied 
e[{N) ^ ai ^ l-A^(A^) 



A'i(A^) - 1 ^ "2 ■ 

It is easy to see that that for the map / to be hyperbolic on N it is necessary 
that X'l > IjAj < 1 and the linearization of / at a;* is hyperbolic with one stable 
and unstable direction. 

Theorem 4.4 Thm. 3] Assume that f is hyperbolic on N . Then 

1. if f''{x) e N for k >0, then lim^^oo f''{x) = x*, 

2. ifvk^N and f{yk-i) = yk for k <0, then limfe^_oo Vk^x^. 

The next theorem shows how we can combine C°- and C^-tools to prove the 
existence of asymptotic orbits with prescribed itinerary. 

Theorem 4.5 J^, Thm. 4] Assume that g is hyperbolic on and f hyperbolic 
on Nq. Let Xg G N„i be a fixed point for g and Xf € Nq be a fixed point for f . 

1- // 

NiJ^ N„, ^ Nm (4.38) 

then there exists xq G Nq such that 

fi-i o fi-2 ° ■ ■ ■ ° faixa) e Ni for i = l,...,m, 
g'' ° fm-i o ■ ■ ■ ° faixa) e N„i for k > 0, 

lim g'' o o • • • o fo{xo) = Xg. 

k—>-oo 
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2. // 

Ni^^ ... Nm (4.39) 

then there exists a sequence {xk)^^_^, f{xk) — a^fc+i for k < such that 

Xk e No for k <0, 

ft-i o fi-2 o ■ ■ ■ ° faixa) ^ N.i for i = l,...,m, 
lim Xk = Xf. 

k — ^ — oo 

3. // 

NiJ^ N„, ^ N„, (4.40) 

then there exists a sequence {xk)1=^oo' fi^k) — Xk+i for k < such that 

Xk <E No, k< 0, 
fi-i o ft-2 o • ■ • o /o(xo) e Ni for i = 1, . . . , m, 
g'^ ° fm-i ° ■ ■ ■ ° fo{xo) Nm for n>0, 

lim Xk = Xf, 

k — > — oo 

lim g'' o /„,_! o • • • o /o(xo) = Xg. 

k — >oo 

The above theorem can be used without any modifications for proving the ex- 
istence of trajectories converging to periodic orbits. In this case we consider 
higher iterates of maps / and g in (4.38), ( 4.39| ) and ( 4.40| ). 



4.2 How to prove an existence of an heteroclinic orbit, 
fuzzy sets. 

To prove an existence of an heteroclinic orbit we want to use the third assertion 



in Theorem 4.5 for g = f, but in order to make the exposition easier to follow 
we use two different maps / and g. Observe that to apply this theorem directly 
one needs to know an exact location of two fixed points a;/ G No and Xg £ Nm, 
because the sets A'o and Nm are centered on Xf and Xg respectively. But exact 
coordinates oi Xf and Xg are usually unknown. We overcome this obstacle in 
three steps as follows 

1. Finding very good estimates for Xf and Xg. In this paper we use an 
argument based on symmetry to obtain tight bounds for Xf and Xg. In a 
rigorous interval Newton algorithm was used. Let us denote by Alf and Mg 
obtained estimates for Xf and Xg, respectively. 

We choose one fixed point Xf £ Mf and Xg G Mg for further considerations. 

2. C^-computations, hyperbolicity We choose a set Uf, Mf C Uf, on 
which we compute rigorously [Df{Uf)]. Then we have to choose a coordinate 
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system, in which the matrix [Df{U /)] will be as close as possible to the diagonal 
one. In this paper we have chosen numerically obtained stable and unstable 
eigenvectors. Let us denote these eigenvectors by u and s, where u corresponds 
to unstable direction and s is pointing in the stable direction. Assume that this 
process gives us a coordinate frame in which 

e[{Uf)e',{Uf) < (1 - mf)){\[iUj) - 1). (4.41) 



From (4.41) it follows easily that there exists ai > 0, Q!2 > such that 

em) (4.42) 



X[{Uf)-l a2 e'^iUf) 

Observe that above inequality speci fies o nly the ratio 0:1/02, hence we can find 
a pair (oi, 02) such that condition ( |4.42| ) and the following condition holds 

Mf + ai ■ [-1, 1] • u + 02 ■ [-1, 1] • s C [// (4.43) 

We now define a h-set A'o by 

Nq = t{xf,aiu,a2s). (4.44) 

Obviously / is hyperbolic on Nq. Observe that the hyperbolicity implies unique- 
ness oi Xf in A'^o- 

We do similar construction for g to obtain = t{xg^ Piu, (32s). 
3. Covering relations for fuzzy h-sets. We have to verify the following 
covering relations 

No^No^ Ni (4.45) 

iV„_l N,n ^ Nra. (4.46) 

As was mentioned above we don't know the h-sets A^o, A^m explicitly, but we 
know that 

A^o e A^o = {t{c,aiu,a2s) \ c € Mf} 

Nm e N„^ = {t{c, PiU, P2S) \ C G Mg}. 

Above equations define a fuzzy h-set, as a collection of h-sets. We can now 
extend the definition of covering relations to fuzzy h-sets as follows. 

Definition 4.6 Let f be a continuous map on the plane. Assume Ni, N2 are 
fuzzy h-sets (collections of h-set) and R is a h-set. 

• we say that Ni M- R iff M R for all M E Ni . 

• we say that R ^ iff R ^ M for all M e Ni. 

• we say that Ni =^ N2 iff Mi ^ M2 for all Mi e Ni and M2 £ ^2- 
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With the above definition is obvious that to prove the covering relations in 
equations (4.45) and ( 4.46| ) it is enough to show that 



Ni (4.47) 
Nm-i ^ Nra. (4.48) 

In practice (in rigorous numerical computations) it is convenient to think 
about a fuzzy h-set N as an parallelogram with thickened edges, hence all tools 
developed to verify covering relations for h-scts can be easily extended to fuzzy 
h-sets. 



5 The Lyapunov orbits. 

In this section we present a computer assisted proof of an existence and hy- 
perbolicity of the Lyapunov orbits around libration points. Hence we realize 



here step 1 and 2 from section 4.2 on our way to the proof of an existence of 
heteroclinic connection. 

As in previous section in the symbol of covering relation we will drop the 

degree part, hence we will use N ^ M instead of ^ M for some nonzero 
w. 

5.1 Existence of Lyapunov orbits. 

Theorem 5.1 Let xi = 0.9208034913207400196, X2 = 1.081929486841799903. 

• There exists a fixed point L\ = (a^i,0) G 0+ for P+, such that 

\x\ - xi\ < r]i = Q ■ IQ-^'^ (5.49) 

• There exists a fixed point Li, = {x2,0) £ 6- for P_ such that 

-a;2| < 7?2 = 10"^^ (5.50) 

Proof: We consider two intervals /i := [xi — rii,xi + rji] x {0} C &+, 
I2 ■= [x2 — ^72,2^2 + V^] X {0} C The location of Xi is schematically shown 
on Figure ^. 

Let us recall, that by TTj, : O — > R we denote the projection onto x coordi- 
nate. With a computer assistance we proved the following 

Lemma 5.2 The maps Pi _|_ : /i — > 6- and Pi _ : I2 — > 6+ are well defined 
and continuous. Moreover we have the following properties 

TTi{Pi^+{xi-m,0)) <0, ^i(Fi,+ (xi+7]i,0)) >0 (5.51) 

TTi{PL^_{x2-m,0)) <0, 7ri(Pl_„(x2 +772,0)) >0 (5.52) 
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Figure 5: The Lyapunov orbits and the location of x*. 



Figures ^ and display rigorous enclosures for Pi (xi ±771 , 0) and Pi _ (x2 ± 
?72,0), respectively. 



Now we are ready to finish the proof of Theorem 5.1. From Lemma p.2| 
and the Darboux property it follows that there exist points x\ G int(/i) and 
X2 G int(/2), such that 



Pi^+{xl,0) = {x\,Q) 
Pi(xl,Q)^{xlQ). 



An application of symmetry properties of Pi/2,± (see eq. (|2.7|)) gives 

P+{xl,G) ^ [x\,Q) 
P_(x^,0) = (a:;,0). 



(5.53) 
(5.54) 



(5.55) 
(5.56) 



5.2 Hyperbolicity in the neighborhood of Lyapunov or- 
bits. 

The goa l of this section is to prove that P is hyperbolic in the sense of Defini- 
tion 4.3 in some neighborhood of points L* and 
Let us define 

ui = (1, 2.5733011), si = (-1, 2.5733011), 
U2 = (1, 2.2817915), S2 = (-1, 2.2817915). 

These vectors appear to be a good approximation for unstable (u^) and stable 
eigenvectors (s^) at L* on the (a;, i)-plane. Observe that R{ui) — ~Si, this is 
in agreement with symmetry of P± stated in equation ( p. 71 ). We will also use 
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Figure 6: Rigorous enclosure of Pi +(xi — ?7i,0) (a box in lower left corner ) 
and Pi _|_ (xi + 771 , 0) (a box in upper right corner) 



(ui, Si) later, as the coordinate directions for good coordinate frame in the proof 
of hyperbolicity of P+ and P_ in the neighborhood of L* . 

Let Hi = t{hi, uj, sj) and Hf = t{hi,u^, s|) for i = 1, 2 denote h-sets on the 
{x,i) plane, where 



hi = (0:1,0), /i2 = (2:2,0) 
ai^S-lQ-i", 02=4-10-10 

si = aisi, si = a2S2 

= 2 • lO^"^?/!, ^ 1.2 • 10-*U2, 
= 2 ■ lO^'^si, s| = 2.8 • 10"'^S2 

We assume that Hl,Hf C e+ and H^,H^ C 6-. Observe that h C Hlc 
Hi and I2 C H2 C where sets li were defined in the proof of Theorem |5.l| . 
Let 

W,: = [-?7^,^7^] X {0}, i = 1,2 (5.58) 



where rji where defined in Theorem 5.1. Let Ui, for i = 1, 2 be given by 



U.,=Hl +W^^{{x+p,i) : (x,i) e ff/,(p,0) e (5.59) 



The choices made in ( ^.57 ) are motivated by the following considerations: 



since we want exploit hyperbolicity of P around L* it is desirable to choose 
stable and unstable directions as Si and Ui. Sets Hf (in fact Ui) will be used 
to establish hyperbolicity around L*, hence it desirable to choose them very 
small, as we need to perform a costly rigorous computation of DP on Ui (a 
C^-computation). Sets Hf are used as a link in the chain of covering relations 
between small Hi and relatively large sets iV^ defined later in Section ^. Since 
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1.046825902266446784 
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Figure 7: Rigorous enclosure of Pi _{x2 — ?72,0) (a box in lower left corner ) 
and Pi _{x2 + ?/2, 0) (a box in upper right corner ) 



the unstable eigenvalue is bigger than 10'^ (see proof of Lemma ^.5[ ), we can 
choose Hf to be of three order of magnitude larger than Hf and still have a 
covering relation between them. 

The following lemma was proved with a computer assistance 

Lemma 5.3 The maps P+ : Ui 0+ and P_ : U2 Q- are well defined. 
Moreover we have 

where intervals Ai, B^, Ci, are given below 

Ai = [695.659,696.1085] Bi = [270.3511,270.4973] 

Ci = [1789.9112,1791.46231] Di = [695.61982,696.12441] 

A2 = [573.3983,573.835] B2 = [251.3098,251.4675] 

C2 = [1308.1679,1309.5201] D2 = [573.3613,573.; 



Using the above lemma and symmetry R we can now prove the following 

Lemma 5.4 There exists exactly one fixed point L\ = (a;J^,0) G Ui for P+. 
Moreover we have \x'{ — xi \ < rji. 

There exists exactly one fixed point L*2 ~ (2^2,0) e U2 for P_. Moreover we 
have — X2\ < ?72- 

Proof: We write down the proof for L*, only. The proof for L2 is analogous. 
An easy computation shows that 

dei{[DP+{Ui)] - Id) <0 



hence from Theorem 4.2 it follows that there exists at most one fixed point for 
P-f in Ui. Since Ii C Ui then we know from Theorem 5.1 that one such fixed 
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point L* = (.i;*,0) £ Ii exists. The estimate for \xl — xi\ was also given in 
Theorem |]l[ I 

Lemma 5.5 There exist R- symmetric h-sets Hi and H2, such that \Hi \ C Ui, 
\H2\ C U2,t L\ G Hi and G ^2 0*^^ ^^le following conditions hold 

1. P_|_ is hyperbolic on \Hi\ 

2. P- is hyperbolic on \H2\ 



Proof: We will proceed as it was outlined in step 2 in section 4.2. First we 



need to fin d a c oordinate frame (via an affine transformation) in which the 



inequality ( ^.4l[) is satisfied for (/ = P+,Uf = Ui) and {f = P_,Uf = 1/2)- 
Form Lemma 5^ it follows that P+ is defined on Ui and P_ is defined on 1/2- 



Observe that the transformation of [DP^{Ui)] {[DP_(JJ2)]) to new coor- 
dinates does not depend on the exact location LI (ij). In new coordinates 
LI = L2 ^ 0, but we have to choose the coordinate directions in Ui and U2- It 
turns out that the vectors (ui,Si) which were used in the definition of H^ are 
good for this purpose, as they are reasonably good approximations of unstable 
and stable directions of corresponding Poincare map. 

A short computation shows that in new coordinates we obtain 

[DP^iUi)] c f^'^' I'A, [DP^{U2)] C (l'^' I'A (5.61) 

\ei,2 '^1,2/ \^2,2 ^2,2/ 

where 

Ai.i = [1391.271,1392.239] Ai,2 = [-0.482,0.485] 

£1,1 = [-0.494,0.472] £1.2 = [-0.483,0.484] 

A2.1 = [1146.751,1147.69] A2,2 = [-0.468,0.47] 

£2,1 = [-0.481,0.457] £2,2 = [-0.468,0.47]. 

It is clear that Ai,2 < 1 < A^^ and e-i, 16^,2 < (1 — Ai,2)(Ai,i — 1). Moreover 
^I'l < 1 < (5.62) 



Ai,i — 1 £1,2 



A2,i — 1 £2,2 
We define Hi for i = 1 , 2 as follows 

Hi=t{L*,aiUi,aiSi), (5.64) 



where ai were given in (5.57). Observe that by construction 

\H,\ C \Hl\ C U,. (5.65) 

This shows that P+ is hyperbolic on |i7i| and P_ is hyperbolic on |ff2|- ■ 
With a computer assistance we proved the following lemma 
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Lemma 5.6 Let Hi and H2 be the h-sets obtained in Lemma 5.1, then 

• ^2 ) ' H2 H2 

Proof: Consider the following fuzzy h-sets 

Hi = t{Wi,aiUi,aiSi), (5.66) 



where Wi are defined by equation ( |5.58| ). We assume that \Hi\ C 6+ and 

1^2! ce_. _ 

Observe that Hi € Hi. The fuzzy sets Hi reflect out lack of knowledge of 
exact coordinates of L* . 



0.9208031 0.9208039 









P+(Hi) 




















Hi ^ 















X 



Figure 8: The Set H^ (the large parallelogram), the fuzzy set Hi (a small set 
in the center) and P+{Hi) (the nearly diagonal segment across Hf) illustrating 

covering relation: Hi Hi Hf. Vertical edges (when in color: red and 
blue) are marked by a bold line. 



The following covering relations were established with a computer assistance 
(see Fig. ^) 

HiJ=^Hi^ Hi (5.67) 
Hi i^H2^H2. (5.68) 



The assertion of the lemma follows now immediately from Def. i.t . 



6 An existence of homo- and heteroclinic con- 
nections for Lyapunov orbits. 



In this section we prove with a computer assistance Theorem 1.1. During the 



proof we define h-sets which will be used later in the construction of symbolic 



dynamics in the proof of Theorem 1.2 
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Figure 9: The location of sets Ni along a heteroclinic orbit in {x, j/)-coordinates. 

6.1 An existence of heteroclinic connection between Lya- 
punov orbits. 

In order to prove an existence of heteroclinic connection between and we 
need to find a chain of covering relations which starts close to L* (begins with 
Hf) and ends close to L2 (with For this sake we choose the sets Ni along a 
numerically constructed, (nonrigorous), heteroclinic orbit in the vicinity of the 
intersection of such orbit with the section 8 (see Fig. ||). Let Ni = t{Xi,Ui, Si) 
be h-sets, where 

Xo = (0.9522928423486199945,1.23-10-^) 

Xi = (0.921005737890425169,0.0005205932817646883714) 

X2 = (0.957916338594066441,0.02191497366476494527) 

X3 = (1.030069865952822683,0.00330658676251664686) 

Xi = (0.967306682018305608, 0.003703230165036550462) 

X5 ^ (1.040628850444842879,0.02317063455298806404) 

Xe = (1.081670357450509545,0.0005918226490172379421) 

Xr = (1.046819673646057103,2.13365065043902489-10-^) 

and 



So = 


(-4-10-^1.45-10-5) 


Uq = 


-i?(so)/10 


Si = 


(-4.5- 10-^1 - 10'^) 


Ui = 


-R{s,)/10 


S2 = 


(-1.2- 10-"^, 2.92- 10-7) 


"2 = 




S3 = 


(-1.05-10-7,2.92-10-7) 


"3 = 


-i?(s3) 


S4 = 


(-1 - 10-7,2.9-10-7) 


Ui = 


-i?(s4)/2 


S5 = 


(-1.44- 10-7,5.8-10-7) 


"5 = 


-i?(s5)/6 


S6 = 


(-1.625 - 10-7,3.75 - 10-7) 


ue = 


-i?(s6)/2 


S7 = 


(-8.3-10-7,2.9-10-6) 


Uj = 


-i?(s7)/5 
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Vectors Si were chosen to be a good approximation of the stable direction at 
Xi. Vectors Ui were chosen as symmetric image of Si, but usually with different 
length. 

Remark 6.1 For our computations to succeed vectors Ui might have been cho- 
sen quite arbitrarily (but not too close to Si). For example, numerical explo- 
rations show that if w e ch oose Ui = aR{si) then we can reduce time of rigorous 
computation (Lemma 6^ ) around 3 times in comparison to Ui parallel to x-axis. 

We assume that iVo, ^2, A^4, A^e C e_ and iVi , iVg , TVg , TVt C 9+. With 
a computer assistance we proved 

Lemma 6.2 The maps 

Pi .+ : UN1UN3UN5UN7 — ^ e_ 

and 

Pi _ : NqU N2U NiU Ne — > 0+ 

are well defined and continuous. Moreover, we have the following chain of cov- 
ering relations 

N, "^^^ Nj ''m^ Hi 

Figure ^ illustrates the chain of covering relations from Lemma |6.2| obtained 

by a nonrigorous procedure. 

Now we are ready to prove the part of Theorem 1.1 concerning an existence 
of heteroclinic connections. 

Theorem 6.3 For FCR3BP with C = 3.03,. ^ = 0.0009537 there exist two 
periodic solutions in the Jupiter region, L\ and L2 , called Lyapunov orbits, and 
there exists heteroclinic connections between them, in both directions. 

Proof: We prove only an existence of the connection from to L2. An 
existence of connection in the opposite direction is obtained by symmetry R. 

From Lemmas 5.6 and 6.2 it follows that there exists the following chain of 
covering relations 

H,^H,^ Hi No N, ^^-^ N, N, N, 

^i^- N, ^=^+ N, N, Hi ^H2^H2 



The assertion follows now from Lemma p. 9 and Theorem 4.5 
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b) 





c) 



d) 





f) 




.0S167I2 1.0468175 
0.1)00593(1 




g) 



h) 



Figure 10: A chain of covering relations, a) H( 



2 



Pl/2.- 



No, b) iVo Ni, c) 



Ni ''^^ iV2, d) N2 ""^^ iVa, e) N3 ""^^ A^4, f) A^4 ""^^ N5, g) N5 ""^^ iVg, 

h) A'^g A^7. These pictures are'nt produced by a rigorous procedure, as we 

checked the covering relations by less direct approach to reduce the computation 
time - see section for details 
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6.2 Homoclinic connection in an exterior region. 



In this section we establish an existence of an orbit homocHnic to (see Fig- 
ure ^l]). For this end we find a chain of covering relations, which starts close 
to L2 passes through the sets located in the exterior region and then ends close 
to L^. For this sake we choose the sets Ei along a numerically constructed, 
(nonrigorous), homoclinic orbit in the vicinity of the intersection of such orbit 
with the section 6 (see Fig. 11). 




Figure 11: Homoclinic orbits to L* and Lyapunov orbits. 



We define the following h-sets Ei = t(Yi,Ui, s^), where 



Yo = 


(-2.08509704964865536, 0), 


Yi = 


(1.160261327316386816, 


-0.1812035059427922688), 


Y2 - 


(1.059527808809695232, 


-0.03871458787165545984), 


Y3 = 


(1.082284499686768768, 


-0.0008090412116073312256), 


Yi = 


(1.046834433386131072, 


-0.00002957990840481726976) , 


Y5 = 


(1.081929798158888576, 


-0.0000007068412578518833152) 



So 


= (-1-10- 


^3- 10-») 


uo = 


-R{so) 


Sl 


= (1-10-^ 


8- 10-») 


Ul = 


-4i?(si) 


S2 


= (-3-10- 


81 -10-8) 


U2 = 


-i?(s2)/10 


S3 


= (-1-10- 


-7,23-10-8) 


"3 = 


-i?(s3)/4 


S4 


= (-1-10- 


35 -10-8) 


W4 = 


-i?(s4)/4 


S5 


= (-1-10- 


'8,22817915- 10-15) 


W5 = 


-i?(s5)/2 



We assume, that Eq,E2,E4 C 8+ and Ei^E^^E^ C 9_. 
With a computer assistance we proved the following 
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Lemma 6.4 The maps 



Pi 

2 

Pi 



,+ 



are well defined and continuous. Moreover, we have the following chain of cov- 
ering relations 



Et 



Pl/2,^ 



El 



E- 



E. 



We are now ready to state the basic theorem in this section. 

Theorem 6.5 For PCR3BP with C = 3.03, /x = 0.0009537 there exists a an 
orbit homoclinic to 



Proof: From Lemmas S.4 and 5.6 it follows that 

^1/2,.- 



^1/2 + ^1/2 - p ^1/2,+ p 



^^5 



(6.69) 



Observe that from the definition of Eq it follows that Eq is i?-symmetric. From 
Corollary 3.14, Lemma 5.6 and equation (6.69) we obtain 



H2 = R{H2) ^ R{Hl] 



^ R[E,) ^'4^- R{E,) 



RiEs) "V^- R{E2) "^^^ R{Ei) '^<^- Eo = R{Eo) 



(6.70) 



From (6.69), ( |6.70 ), Lemmas ^.6| and Theorem 4.5 we obtain an orbit homoclinic 
to iS. I 



6.3 Homoclinic connection in interior region. 

In this section we establish an existence of an orbit homoclinic to (see Fig- 
ure |l^). For this sake we find a chain of covering relations, which starts close 
to LI passes through the sets located in the interior region and then ends close 
to LJ. For this sake we choose the sets Fi along a numerically constructed, 
(nonrigorous) , homoclinic orbit in the vicinity of the intersection of such orbit 
with the section Q (see Fig. 11). 

We define the following h-sets Fi — t{Zi, Ui,Si) where 

Zo = (-0.6160415155975000064,0), 

Zi = (0.84668503722876047360,0.17563753764246766080), 

Z2 = (0.94793695784874987520,0.01522141990729746432), 

Z3 = (0.92067611200358768640,0.00032764933375860776), 

Zi = (0.95228425894935162880,0.00001048139819208300) 
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and 



So 


= (- 


1 • 10-^25-10-8) 


uo = 


-R(so) 


Si 


= (1 


■ 10-^92 • 10-9) 


Ul = 


-2.2R{si) 


S2 


= (- 


25. 10-9, f -10-8) 


U2 = 


-i?(s2)/5 


S3 


= (- 


1 • 10-^26-10-8) 


U3 = 


-i?(s3)/6 


S4 


= (- 


1 • 10-^37-10-8) 


U4 = 


-i?(s4)/6 



We assume, that Fq, F2, F4 c 0_ and Fi, F3 c 6+. 
With a computer assistance we proved the following 



Lemma 6.6 The maps 



Pi_ : FoUF2UF4-^e+ 
Pi , : FiUF^^O- 



are well defined and continuous. Moreover, we have the following covering re- 
lations 



Fa "^^^ F4 



Hi 



We are now ready to state the basic theorem in this section. 

Theorem 6.7 For PCR3BP with C = 3.03, /i 0.0009537 there exists a an 
orbit homoclinic to L\. 



Proof: From Lemmas 6.6 and 5.6 it follows that 



Fo Fi F2 F3 F4 Hi ^ Hi (6.71) 

Observe that from the definition of Fq it follows that Fq is i?-symmetric. From 



Pl/2 



1/2; + 



Corollary 3.14, Lemma 5.6 and equation (6.71) we obtain 

Pl/2 



Hi = R{Hi) 4^ R{Hl) i?(F4) i?(F 



3; R{F2) "'4^ 



R[Fi) "V^^ i?(Fo) = Fo (6.72) 



From (6.71), (6.72), Lemmas 5.6 and Theorem 4.5 we obtain an orbit homoclinic 
to LI. I 



Proof of Theorem 1.1: We combine together Theorems 6.3, 3.5 and 6.7 
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7 Symbolic dynamics on four symbols 



The goal of this section is to give a precise meaning and a proof to Theorem 1.2 



As in previous sections in the symbol of covering relation we will drop the 

degree part, hence we will use N ^ M instead of ^ M for some nonzero 
w. 



From Lemmas 5^ and 3^ we know that there exists the following chain of 
covering relations 



Pl/2,- 



No 



Pi/ 



Pl/2,- 



^1/2,4 



Pl/2.- 



Pl/2A 



Pi/ 



Ho 



Ho 



From Lemmas 
From Lemmas 



5J and ^T2| we have R{H{) = ff, for i = 1, 2. 
5.6, 6.2 and Corollary |3.14 it follows that 



H2 = R{H2) ^ R{Hl 



R{N^) R{Ne) ''<^- R{N5) 

"4^^ RiN,) "V^- RiN,) RiN,) ^'4^- R{N,) 



RiNo) 



RiHf) ^ R{Hi) = Hi 



From Lemma 6.4 and the proof of Theorem |6.5| it follows that 

Pl/2.- ^ Pl/2,i 



^1/2 + p Pl/2 - ^1/2 + p 



Ho 



(7.73) 



(7.74) 



(7.75) 



and 



H2 — R(H2 



RiHl) ^ R{E,) 



R{Ei 



1/2,+ 



RiEs) "^'4^- R{E2) "4^+ R{Ei) ''k^- Eo = RiEo). 



From Lemma |6.6| and the proof of Theorem 6.7 it follows that 



(7.76) 



Er 



Pl/2,- 



Ei 



Pi/ 2;+ „ Pi/ 2.- 

=^ r2 



E. 



Pl/2,^ 



E4 



Pl/2,- 



Ht 



H, 



(7.77) 



and 



Hi = R{Hi) ^ RiHl)''^^^ 



R{Ei 



Pl/2.- 



RiEs) 
RiEi 



RiE2) 

"i?(Fo) = Fo. (7.78) 



We will construct now the symbolic dynamics on four symbols. The con- 
struction is a little bit involved, because we have four different maps in all 
covering relations listed above. 

We assign symbols as follows: 1 - the set Hi, 2 - H2, S - Eq and 4 - Eq. The 
covering relations allow for transitions 1 ^ 1, 1 — > 2, 1 — > 4, 2 — > 1, 2 — > 2, 
2 ^ 3, 3 ^ 2 and 4^1. For each such transition i ^ j we associate a pair 
This defines a set of admissible pairs F. 
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For any (a, /?) G F we define a map f{a,i3) as follows 



/(a,/3) 









if (a,/3) = 


(1,1), 


f"- o ^"1/2,+ 


iPl/2,- 


0Pl/2,+ )^0P+ 


if (a, /3) = 


(2,1), 


^+ ^1/2 - 


iPl/2,+ 


A/2,-)*oP_ 


if {a, (3) = 


(1,2), 








if 


(2,2) 


' ° (^1/2,- 


- ° ^1/2, 


+ )'°A/2,+ 


if (a,/3) = 


(4,1) 


P^,- {Pl/2.- 


^- ° A/2, 


_)2oP2 


if (a, /3) = 


(3,2) 


P+ Pl/2,- 


(A/2,+ 


Pl/2,0)' 


if (a,/3) = 


(1,4) 


.^"^oPl/2,+ 


(A/2,- 


oPl/2,+ ? 


if {a,f3) = 


(2,3) 


{1,2,3,4}^ be 


defined 


as follows c G 


Sr iff for 


every i 



id, Ci+i) e F. We set Si = Hi, S2 = H2, S3 = Eq and S4 = Fq. 

We can now formulate the theorem about an existence of symbolic dynamics 
on four symbols symbols 

Theorem 7.1 For any sequence a = {ui} € Sr there exists a point xg £ S^o, 
such that 

• its trajectory exists for t G {—00, 00) 

• Xn= /(a„,a„_i) O ■ • • O f(a2,at) ° /(ai,ao)(^o) G Sa„ for n > 



Moreover, we have 



-1 



-1 



(0-1,0-2) •'(ao,a-l) 



Jxq) G Sa„ for n<0. 



periodic orbits: // a is periodic with the principal period equal to k 
can be chosen so that Xk = Xq, hence its trajectory is periodic. 



th- 



en Xo 



homo- and heterocUnic orbits: // a is such that ak = i- 
Qffc = i-i- for k > k^, where i+ G {1, 2}, then 



for k < k- and 



lim Xn = L* 



n — ^ — oc 



lim Xn — L] 



Proof : From chains of covering relations (|7J3|) , (|7J^) , ( [ZJSI) , ([tJGI) , (fzJTl) , 
( 7.78 ) and Theorem 3^ we obtain the statement on periodic points for periodic 
a. To treat a nonperiodic a we approximate it with periodic sequences /3„ with 
increasing periods to obtain sequence of points and after eventually passing 
to a subsequence we obtain xq with desired properties. 

The statement on homo- and heteroclinic orbits is an easy consequence of 
Theorem 4.5 and the hyperbolicity of P± on Hi established in Lemma 5.5 I 

Our methods do not allow to make any claims about the uniqueness of xo 
for a given a. The only claims of this type we can make is if an = « for al\n £ "Z 
then Xo = L*. 
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8 Numerical aspects of the proof 



In this section we give details of the computer assisted proofs of Lemmas 5.2 



5.3, 5.6 and 5.2. As in previous section in the symbol of covering relation we 



will drop the degree part, hence we will use N =^ M instead of N M for 
some nonzero w. 

8.1 The existence and continuity of Poincare maps. Hy- 
perbolicity on Uj. 

All proofs required to check first that suitable Poincare maps {P±, Pi ±) are 
defined on some parallelograms (supports of our h-sets) on Q±. For this end the 
parallelogram, Z , was represented as a finite union of small parallelograms, Zi, 
and each of Zi^s was used as an initial condition for our routine computing the 
necessary Poincare map, Pi/2.± or P±. We divided horizontal edges on n equal 
parts (a horizontal grid) and vertical edges on m equal parts (a vertical grid) 
and hence we covered Z hy n x m parallelograms. Our routine was constructed 
so that, if completed successfully, then we can claim that Zi is contained in the 
domain of P and the computed image contains P{Zi). Our routine is based on 
the C° and C^-Lohner algorithms @ |3). 
We had to prove the following assertions 



1. (in Lemma x2) Pi ^ is well defined and continuous on Ii and Pi _ is well 



defined and continuous on l2- 



(in Lemma 5.3) P+ is well defined and smooth on Ui, P- is well defined 
and smooth on C/2. 



(in Lemma |5.6| - equations ( ^.67| ) and ( 5.68| )). P_ is well defined and 



continuous on Observe that since Hi C C/i, H2 C U2, then the 

previous assertion guarantees an existence and continuity of P+ on Hi 
and P_ on H2- 



4. (in Lemmas 3.2, 3.4 and |6.6| ) Pi is well defined and continuous on 
Ni, N3, N5, Nt, Eq, E2, Ei, Fi, F3 . Pi _ is well defined and continuous 
on Nq, N2, N4, Nq, El, E3, Fq, F2 and P4. P_ is well defined on E5. 

The first assertion follows easily from the second one. We reason as follows: 
since Ii C Ui, then an existence of P_ (P+) on Ii {I2) implies that also Pi _ 
(P^ ^) is defined. 

To prove the second assertion we cover Ui by finite number (13 x 13) of 
parallelograms. Then we compute an image of each part and an enclosure of 
the derivative of the Poincare map using a routine based on C^-Lohner algorithm 
recently proposed in |23]. As a consequence we obtain an estimation of DP± 



(see Lemma 5.3). Parameter settings used in these computations are listed in 
Table |^. Let us stress also, that a successful termination of our routine proves 
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set 


order 


step 


horizontal grid 


vertical grid 


Ui 


5 


0.007 


13 


13 


U2 


5 


0.007 


13 


13 



Table 1: Parameter settings of the Taylor method used in C^-computations - in 



the proof of Lemma 5.3 



also that P+ and P_ are defined on Ui and U2, respectively. From the standard 
theory it follows that P± are smooth on their domain. 

To prove the third and fourth assertion we proceed in the similar way. We 
cover each set by finite number of parallelograms and compute an image of each 
parallelogram. Since an estimation of the derivative of the Poincare map is not 
necessary we have used a C°-Lohner algorithm |lj, . Parameter settings for 
these computations are listed in Table H. 



covering relations 


order 


step 


h. grid 


V. grid 


H, ^ Hi "^i^^ No N, ''^^ N2 

N2 N, N, N, 
N, ''^^ N, ^i^- N, HI ^ H2 , 


5 


0.01 


1 


1 


fl/2 + p Pl/2- ^1/2 + p 


5 


0.02 


1 


1 


^0 =^ ^1 =^ -f^2 -Ta 
""^^ F, ""m-- Hi ^ H, 


5 


0.02 


1 


1 



Table 2: Parameters of the Taylor method used in the proof of an existence of 



the Poincare map in Lemma p. 2, Lemma 5.6, Lemma 5.4 and Lemma 13.6 



8.2 Details of the proof of Lemma |5.2 



To prove inequalities (5.51),(53^) we had to compute rigorous enclosures for 



J- i + ± 771, 0) and Pi _{x2 ± 772, 0), respectively. The values of the time step 
and the order of the Taylor method used in our routine are listed in Table ^. 
Figures ^ and ^ display the actual enclosures obtained. 
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orbit 


order of Taylor method 


time step 


LI 


20 


0.05 


LI 


19 


0.055 



Table 3: Settings used in the proof of (p.5l|) (L^f-row) and (5.52)(L2-row) 



8.3 How do we verify covering relations - details of proofs 
of Lemmas |6.2| , |6.4| and |6.6 



This is the most computationally demanding part of our programm. 

In principle the same rigorous computations can be used to obtain both an 
existence of Poincare maps and covering relations, but in practice this does'nt 
work, i.e. it will result in an enormous computation time (see the discussion in 
Sec. 6 of 1)). 

It turns out that once an existence of Poincare map is established, we can 
reduce the computations to the boundary of our h-sets and one interval inside. 



only (see Lemmas 3.3 and 8.4). Now, when we compute an image of an edge 
/, we still have to divide it into subintervals, but the number of subintervals of 
the order of square root of the number parallelograms need to achieve the same 
accuracy on the parallelogram build on two intervals of the linear size similar 
to that of /. 

In order to establish an existence of covering relations we need to verify the 



assumptions of Theorem 3.9 



To facilitate a discussion about various conditions implying Theorem 3.6 we 
introduce the following 

Definition 8.1 Let / ; — > he a continuous map and let Ni = t{ci, ui, si) 
and N2 — t{c2,U2, S2) be two h-sets. 

We say that f satisfies condition ah, aO, a, b+, b- on Ni and N2 if 

ah: there exists qq G [—1,1] , such that 

/(cat, ([-1,1] X {qo})) C int {Ni U ^UN^) 

aO: /(|iVi|) niV2+ = 

a: /(|7Vi|) Cint(iV^U|iV2|U7VJ) 

b+: fiNi") C int(Ar^) and f{N{^) C int(A^2'') 

b-: /(iVl'^) C int(iV2'') and fiN[^) C int(Ar^) 

We say that f satisfies condition b on Ni and N2 if either b+ or b- is 
satisfied. 
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Remark 8.2 

Theorem 



Obse 



that conditions (3.2c ) 



(P.24 ), ^3.2di ) and l[3.26J from 
3.!\ coincide with conditions ah, aO, b+ and h-, respectively. 



Observe that condition a implies conditions ah and aO. 

The following lemma gives sufficient conditions for an existence of covering 
relations for injective maps. 

Lemma 8.3 Let / : — )■ be a continuous map and let Ni and N2 be two 
h-sets. Assume that f is injective on \Ni\ and f satisfies condition b on Ni 
and N2 and the following condition a' 



f{d\Ni\) C int(iV^ U \N2\ U N^). 



The 



Proof: From Remark 8.2 and Theorem 3.9 it follows that it is enough to verify 
condition a. This follows easily from a' and the Jordan theorem (see ||], page 
180). I 



Figures illustrating covering relations obtained in Lemmas 3.2 and 5.6 sug- 
gest that condition a is satisfied in all relations. Unfortunately the verification 

of condition a ( or a') pose the following difficulty: In the relation iVi A'2 
the set |iVi| is mapped across of without touching its horizontal edges, but 
if |A^2| is small then we need a very good estimation of image of horizontal edges 
of Ni . This forces us to make a very fine partition of the boundary of A^i , take 
small time steps and a high order in the numerical method resulting in very 
long computation times. 

Above phenomenon is illustrated on Fig. which shows enclosures obtained 
from our rigorous routines. On this picture we can see rigorous enclosure for 
an image of Pi/2,-{dNQ). This image was obtained as follows: we divided the 
boundary of the set Nq into some number of subintervals (see Table ^ and 
computed an image of each part via Pi/2,-- This picture shows, that much 
tighter enclosure of an image of horizontal edge was required compared to an 
enclosure for an image of vertical edges (for example edge {NqY^ was divided 
into 8 equal parts, but (Nq)^^ into 5 equal parts). In other covering relations 
this disproportion was often much bigger. 

To deal with this problem we use the following lemma, in which we indirectly 
verify conditions aO and ah instead of a' . This approach allowed us to reduce 
the computation time by a factor of 5. 

Lemma 8.4 Let Ni = t{ci,ui, si), N2 = t{c2,U2,S2) are h-sets, f : |iVi| — > 
an injection of class C^. Let j be a horizontal line in \Ni\ connecting vertical 
edges given by 



7: [-1,1] 3t — >ci+t-uie \Ni\ 



.79) 



Let g — (31,52) = 0/07. Assume f satisfies condition b on Ni and N2 
and the following conditions hold: 
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al ^(t)^0/orte(-l,l) 

a2 there exists to € (—1,1) such that /(7(fo)) S int(|Ar2|) 
a3 f-^N+)n\Ni\ = ^ 

Then Nid^ N2. 

Proof: We need to show that conditions ah and aO are satisfied. 

Observe that condition aO follows immediately from condition a3 and in- 

jectivity of the map /. Namely, by applying / to both sides of a3 we obtain 

^2+n/(|iVi|) = 0. 

We now show that condition ah is true. 

For this end wc consider / o 7 in the coordinates induced by the map Cjva- 
In these coordinates 

I7V2I = [-1,1] X [-1,1] (8.80) 

= [1,00) X (-(X),(X)) (8.81) 

= (-00, -1] X (-CX), 00) (8.82) 

/07 = g (8.83) 

Without any loss of generality we can assume that 

^(i)>0, for te{-l,l), (8.84) 

Hence gi is a strictly increasing function and from condition b it follows that 
b+ is satisfied. 

We define two numbers 

t* = min{t > to I fhit)) e d\N2\}, (8.85) 

t, = max{t < to I /(7(t)) e d\N2\}. (8.86) 

From conditions a2, aS and b it follows that these numbers are well defined 
tt, < t* and 

/(7([-l,i*)))cint(iV^), (8.87) 

/(7((r,l]))cint(iV20, (8.88) 

/(7((i*,i*))) Cint(l7V2|). (8.89) 

To finish the proof observe that from condition b it follows that 

/(7(±l))Gint (N^UN^) 



36 



Remark 8.5 Observe that Cj^^{x) — A ^{x ~ C2), where A = [uj,s|']. 
Hence 

^(t) = ^ Ar/d/(7W).,ui,„ (8.90) 

where ui = (ui.i,ui^2)- 



Remark 8.6 Observe that, when Ni — t{ci, u, s) and N2 = t{c2, au, S2), which 
means that the unstable coordinate direction for both h-sets coincide, then 



dg 



dt 



'-{t)^a{A-' ■df{j{t))-A), 



.91) 



Hence it is enough to look at the (1,1) entry of df expressed in the coordinates 
of the h-set N2. 



1.046816 



Figure 12: An example of the rigorous enclosure of an image of ONq in relation 



N7. 



In Table ^ we present settings used in the proof of Lemma 6.2. In particular 
the parameter grid gives the number of equal intervals into which we divide an 
edge. A positive time step means that Lemma 3.3 is used to verify a covering 
relation. A negative time step means that we use Lemma 3.4 and symbolizes 



the fact that we compute an inverse of the Poincare map to verify condition 
a3. Parameter settings for the verification of al and a2 are given in Table ^. 
In this table order (m) and step(m) denote an order and a time step of Taylor 
method which we use to prove al and order (c), step(c) denote an order and 
a time step of Taylor method which we use to prove a2. Parameter grid(m) 
denotes a number of equal intervals used to cover the curve 7 in condition al. 
To verify condition a2 we usually compute image of the center of the set 

{to = 0). Only in the case N7 ^=>^ i7| we used to = 0.228. 



37 



covering relation 


edges 


grid 


order 


step 




(H^)''^ and (i7i2)te 


6 


4 


0.01 


m ''^^ No 


and {HfY^ 


7 


5 


0.01 




N^" and TVf 


1 


5 


-0.01 


No ''i^" iVi 


iVo'''^ and Nl," 


40 


8 


0.04 




Nl" and TV^'^ 


8 


6 


-0.01 


N, "^^^ N, 


Nl" and iVl'^ 


25 


5 


0.01 




iVl'^ and TVf 


6 


6 


-0.004 


N2 N3 


iVa'''^ and iV^"^ 


5 


6 


0.004 




iVl'^ and TVf 


3 


6 


-0.005 


Ns "^^^ N, 


iVg'''^ and iYl"^ 


5 


6 


0.004 




Nl" and TVf 


2 


6 


-0.01 


^4 ""^^ N5 


iVX'^ and Ni^ 


15 


6 


0.006 




Nl^ and TVf 


2 


6 


-0.01 


N, N, 


iVg'-'^ and Ni"" 


32 


6 


0.01 




Nl" and TV^'' 


8 


6 


0.01 


Ne ''^^ Nj 


and iV^'^ 


5 


6 


0.01 




(iJl)'"^ and {Hiy" 


2 


5 


-0.01 


N, Hi 


N^" and iV^"^ 


33 


5 


0.01 



Table 4: Parameters of the Taylor method used in the proof of Lemma 3.2 



covering relations 


order(m) 


step(m) 


grid(m) 


order(c) 


step(c) 


No ""m^- N, ""^^ N, 


4 


0.01 


1 


6 


0.01 


N2 ""^^ iVa 


4 


0.004 


4 


6 


0.004 


7V3 ""^^ iV4 


4 


0.003 


2 


6 


0.003 


7V4 ""^^ iVs 


4 


0.004 


1 


6 


0.004 


iVs ''^^ No "^^^ Hi 


4 


0.01 


1 


6 


0.01 


HI ^ H2 


4 


0.01 


1 







Table 5: Parameters of the Taylor method used in the proof of conditions al, 
a2 in Lemmas 3.2 and 5.6 
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covering relation 


edges 


grid 


order 


step 


Eo ''^^ El 


E'l" and El" 


33 


7 


-0.04 


El" and E'q" 


12 


7 


0.04 


Pl/2 - 

El =4- E2 


El" and 


12 


8 


-0.02 


El" and 


25 


7 


0.02 


^1/2 + 


and El" 


1 


7 


-0.02 


E^" and E;^'^ 


50 


7 


0.03 


-P]^ / 2 . — 


and El" 


2 


7 


-0.02 


El" and i^l'^ 


7 


7 


0.03 


i?4 ''^^ E, 


El" and i;^'^ 


2 


7 


-0.02 


and E't 


20 


7 


0.03 


E5 ==^ H2 


(Hir and 


1 


5 


-0.02 


£;,^'= and El" 


4 


5 


0.01 



Table 6: Parameters of the Taylor method used in the proof of covering relations 
in Lemma 6.4. 



covering relations 


o(m) 


s(m) 


g(m) 


0(C) 


s(c) 


J'1/2 + ^1/2 - p ^1/2,,+ p 


4 


0.02 


1 


7 


0.02 




5 


0.02 


1 


8 


0.02 


Fo Fi 


5 


0.01 


2 


8 


0.02 


Fi ''^^ F2 F, ''^^ F, Hi 


4 


0.01 


1 


8 


0.02 


Hi ^ Hi 


4 


0.01 


1 







Table 7: Parameters of the Taylor method used in the proof of conditions al, 



a2 in Lemma 3.4 and Lemma 
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Parameters of the Taylor method used in the proof of conditions of Lemma 6.4 

are hsted in Tables || and 0. In the relation ^=4> if| we do not compute an 
image of the center of set but an image of the point C — — 0.155u5. 

Parameters settings used in the proof of Lemma 6.6 are listed in Tables |^ 
and 0. 



covering relation 


edges 


grid 


order 


step 




F^" and Ff 


50 


6 


-0.01 


Pl/2 - 

Fo ^ F, 


FS" and F^'^ 


20 


8 


0.02 




F^" and F^^ 


9 


7 


-0.02 


F\ ==^> F2 


Fl^ and 


330 


7 


0.02 




F^" and Fg*^ 


1 


7 


-0.02 


Pl/2.- 

F2 ^3 


F^" and F!f 


35 


7 


0.03 




Fl" and Ff 


1 


7 


-0.02 


F3 ''^^ F, 


F^^ and Fl'' 


10 


7 


0.03 




(Hi)'"' and 


3 


7 


-0.02 


F, Hi 


F^" and Fi" 


45 


7 


0.03 




Hi" and iJf 


3 


8 


-0.02 


Hi ^ H, 


and [Hiy^ 


7 


7 


0.03 



Table 8: Parameters of the Taylor method used in the proof of covering relations 
for sets Fi. 



8.4 Verification of covering relations for fuzzy set - details 



of the proof of Lemma 5.6 



In this subsection we discuss how we verify covering relations for fuzzy h-sets. It 
is convenient to think about a fuzzy h-set N as an parallelogram with thickened 
edges. Wejiefine the support, left and right edges and left and right sides of a 
fuzzy set N as follows 

where by (Z) we denoted a convex hull of the set Z . We introduce one more 
notation for allowed image of the h-set covering N 



strip(7V) = Pi int (M' U \M\ U W) 



40 



Lemmas B.3 and 3.4 can be easily adopted to fuzzy h-sets. Namely we have 

Lemma 8.7 Let / : M'^ ^ be a continuous map and let Ni and N2 be two 
fuzzy h-sets. Assume that f is infective on \Ni\ and the following conditions af 
and bf are satisfied: 

af /(l^il) Cstrip(7V2) 

bf either f{N[^) C int(7V^) and f{Nl^) C uii{N^) 
or /(iV{^) C int(7V2'') and fiN^") C int(7V^). 



Then Ni 



Lemma 8.8 Let Ni = t{Wi,ui, si), N2 = t{W2,U2,S2) be fuzzy h-sets and 
f : [7Vi| — > an injection of class . Let j be a fuzzy horizontal line in \Ni\, 
given by 

7 : [-1,1] X 1^1 3 {t,wi) — >wi+t-uie \Ni\ (8.92) 
For W2 S W2 let N2^W2 — t{w2,U2, S2) and 

9w2it,Wl) = {gw2,l,9w2,2){t,Wi) = c~^l^^ o f Oj{t,Wi). 

Assume that Wi is connected, f satisfies condition bf on Ni and N2 and 
the following conditions hold: 

afl wi) ^0 forte (-1, 1), wi G Wi, W2 £ W2 

a£2 there exists to G (—1, 1) and Wi G Wi such that /(7(ioi ^^^i)) G int(|A^2|) 

af3 f-\N+)n\N,\^(D 



Remark 8.9 Observe that (compare Remark 8.1) that ' ^ does not depend 
on W2 and is given by formula ( 8.9(\ ) with the matrix A depending only on U2 
and 32- 

Let us describe how above conditions afl, af2, af3 were verified for relations 
( ^.67 ) and (5.68), which we rewrite below for the convenience of the reader 



Hi 



Hi 
H2 



■Ht 
H2. 



(8.93) 
(8.94) 
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Let us recall (see Section 5.2) that the h-sets entering above covering relations 
are given by 

Hi = t(Wi,aiUi,aiUi) 

Hf = t(/ii,2-10-^ui,2-10-^si) 

Hi = t(/i2,1.2- 10"%,2.8- 10~*S2) 

For all covering relations, N Af, listed above the unstable vectors for 
both h-sets entering the relations are proportional, hence we can apply Re- 



mark 8.6 and look on (1, l)-entry of dP± expressed in CM-coordinates. 



Observe that, since Hi C Ui then from Lemma 5.3 and equation ( ^.61 ) we 



obtain an enclosure for [DP±{\Hi)\ expressed in iJi-coordinates. From inspec- 
tion of Ai l it follows that condition afl holds for covering relations Hi Hi 

and Hi =^ H^. Since is not contained in U2 we had to verify condition 

P_ ^ 

afl for the relation iJ| => H2- Parameter settings for this computation is 
added as the last row to Table ||. 

Condition af2 is c learl y satisfied with W2 ~ L\ for relations ( ^.67|) and 
1^2 = L2 for relations (5.1 



Conditions af3 and bf must be verified in direct computations. Parameter 
settings for these computations are given in Table ^. 

It turns out that some inclusions involved in condition b can be verified at 
the same time. For example, to prove that Hi Hi =^ Hf we need to verify 
condition b+ for both relations. 

P+iHDcHl, P+iHi'^)cHi, (8.95) 
P+iEDciiHly, P+{H[^)cz{Hiy. (8.96) 

It is sufficient to show that 

P+iHl") C HI n {Hfy (8.97) 
P+{H'^) C H'lDiHfy (8.98) 

P_ — p_ 

Similarly, to prove that i/| H2 => H2 we must verify condition af3 for 
both relations. Namely, we have to check that 

P^\H+)n\H2\=(ll, (8.99) 
Pl\H+)n\Hl\^9. (8.100) 



Since |iJ2| C (compare ( 5.57 and ( 5.66 )) it is sufficient to show ( S.lOO ) 

only. 

9 The pcr3bp program. 

The program pcrSbp is based on CAPD C++ package Q developed at the 
Jagiellonian University, Krakow, Poland. A whole proof took less 32 minutes 
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covering relations 


edges 


grid 


order 


time step 




(^1)'"= and {HiY" 


2 


6 


-0.01 




(HiY" and {HiY" 


3 


6 


0.01 


^ Hi 


(Hi)'"' and 


2 


5 


-0.01 




02)''^ and (i?2)*" 


4 


8 


-0.02 


Hi =^ H2 


{HlY'^ and 


32 


5 


0.01 


H2 H2 


and (^2)'" 


2 


8 


0.02 



Table 9: Parameters of the Taylor method used in the proof of the covering 
relations for fuzzy sets. 



on Celeron l.lGhz processor (compiler Borland C++ 5.02 was used in these 
computations). The source of our program can be found at iQ. Below we 
briefly describe basic data structures used by our program. 
We used the following basic classes: 

class interval - realizes an interval arithmetic (see [|l^ and references given 
there). An interval is a a pair of real numbers, each of the double real type 
(64 bits). 

class intervalVector - realizes interval arithmetic on vectors 

class intervalMatrix - realizes interval arithmetic on matrices 

class set - an abstract class to handle various representations of subsets of M" 
used in C°-Lohner algorithm @, 0, ||. 

class rec2set - a class derived from the class set. This is a class of subsets in 
M" of form x = Xq + Br + Crg (doubletons in the terminology introduced 
in0). 

class clset - an abstract class used in C^-computations via C^-Lohner algo- 
rithm ||. 

class clrect - a class derived from the class clset. This class represents subsets 
of R" a,s X = xq + Br and the derivative of the flow with respect to initial 
conditions as D = Dq + BqR. 

class clrect2 - a class derived from the class clset. This class represents 
subsets of M" a,s x = xq + Cro + Br and the derivative of the flow with 
respect to initial conditions as 13 = Dq + BqR. 

class hSet - an abstract class representing h-set in M". 
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class TripleSet - a class derived from the class hSet. This class implements 
a special case of the hset, an h-set on the plane with one stable and 
unstable directions. This class is used to automatize a verification of 
covering relations. 

class Taylor - a class implementing one step of a both C"- and C^-Lohner 
algorithms for ODEs 23 1. It provides the tools which allow to com- 



pute rigorously the evolution of sets described by the classes given above 
along the trajectory of an ODE. The Taylor coefficients of the solutions 
of an ODE are generated using an automatic differentiation (see and 
references given there). 

class PoincareMap - based on the class Taylor provides the tools to compute 
the Poincare map and its Jacobian matrix for an ODE for set classes 
described above. 

The basic methods and operators of the class PoincareMap are: 

set_order(int) - sets an order of the Taylor method in the Lohncr algorithm. 

set_step(interval) - sets a time step of the Taylor method in the Lohner al- 
gorithm. 

intervalVector operator () (set *) computes an image of an abstract set via 
Poincare map. 

intervalVector operator ()(clset *, intervalMatrix &der) - computes an 
image of an abstract clset. On output in der parameter contains a matrix 
of partial derivatives with respect to initial conditions of the fiow, ^ {t, x) 
for all X defined by parameter clset and for t G [ii,t2], where [^1,^2] 
contains Poincare return times for all initial conditions defined by clset. 
From this matrix we can later extract easily ^ (see details). 

Below we provide a guide to the main functions in our programm. The main 
function is located in the file pcrSbp.cpp. Other important functions are: 

Ll_periodic_orbit() and L2_periodic_orbit() located in the file lyapunov.cpp. 



These functions prove the inequalities in the Lemma ^.2| , from which an 
existence of the Lyapunov orbits follows. 

domain() located in the file domain, cpp, used to verify an existence of the 



Poincare map for the chain of covering relations form Lemma 6.2 



domain_Ll() and domain_L2() located in the file domain. cpp, used to verify 
an existence of the Poincare map for covering relations Hf =^ Nq and 



iJ| H2 respectively. ( in Lemmas |6.2| and 5.2) 



domain Jnter lor and domain_exterior() located in the file domain. cpp, 
used to verify an existence of the Poincare map for covering relations 
Fo =^ . . . =^ F4 =^ H"} =^ Hi and Eq =^ . . . =^ E5 =^ i?| (in 



Lemmas 6.6 an 6.4) 
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hyperbolic_Hl() and hyperbolic_H2() loc&ted in file hyperbolic. cpp. These 
functions computes DP on the sets Hi and H2 (see Lemma 5.3) and proves 
covering relations Hi Hi and ff| =^ H2 from Lemma 5.6. 



cover _heteroclinic() located in the file cover, cpp, contains a proof of Lemma 6.2 
namely it verifies the relations H^ A^o =^ • • • =^ ^7 



Hi 



cover_interior() and cover _exterior() located in the file homoclin.cpp, con- 
tains a proof of covering relations from Lemma 6.t and Lemma 6.4. 



monotone(), monotone Jnterior() and monotone_exterior() - located in 
the file monotone. cpp. These functions performs a verification of the 
monotonicity on the horizontal central lines connecting vertical edges in 
h-set. These computations are needed when we use backward computa- 



tions. This happens in Lemma B.4 and Lemma B.8 



10 Concluding remarks, future work 

There are several directions in which this research can be extended. 

First, all the methods presented in this paper are not restricted to the par- 
ticular parameters of the Oterma comet, other parameters may require slight 
changes in the definition of the sets on which covering relations are build, but 
the method itself will be the same. Basically this method can be applied to 
prove a symbolic dynamics in any system for which numerical simulations indi- 
cate an existence of some kind of hyperbolic behavior, for example here we have 
homo- and heteroclinic chains. 

Another interesting problem is the question of an existence of a hyperbolic in- 
variant set claimed in , where the authors assumed an existence of transversal 
homo- and heteroclinic connections between Lyapunov orbits and then followed 
the standard dynamical system theory argument from the Birkhoff-Smale homo- 
clinic theorem. Since we did'nt computed here unstable and stable manifolds, 
we cannot use these arguments. Observe also that an rigorous computation 
of stable and unstable manifolds for our problem appears to be very difficult 
(requires very extensive C^-computations). Hence developing tools which avoid 
a direct computation of invariant manifolds is of interest. In this context we 
formulate the following conjecture. 

Conjecture 10.1 Let f be a diffeomorphism. Let Nq, Ni be h-sets. Assume 



that f is hyperbolic on Nq and Ni (in the sense of Def.4_^). Assume that we 
have the following sequences of covering relations 

Nod^Nod^Aid^A2^---^As^Ni 
Ni Ni Bi B2 =U ■ ■ ■ ^ Br ^ No, 
then there exists k >l and S C jA^ol U \Ni\, such that 



f'^iS) = S, i.e. S is an invariant set for f^ 
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• S is hyperbolic (in the standard sense - see for example ^) 

• the map tt : S* ^ E2 = {0,1} given by TT{x)i = j iff /'^^{x) G \Nj\ is 
one-to-one. 



Another interesting problem is a question of stability of obtained results with 
respect to various extension of PCR3BP. By this we mean the following: 

• Does the symbolic dynamics persists if the Jupiter orbit become an ellipse 
with a small eccentricity (which is the case in nature)? This can be seen as 
a small periodic perturbation to the ODE describing PCR3BP. We believe 
that an answer is positive. Obviously in this context one can consider a 
more general question: 

Assume that we obtained a symbolic dynamics for an ODE x' f{x) us- 
ing covering relations. Does this symbolic dynamics persists for an nonau- 
tonomous ODE x' = f{x) -\-e{t,x) ife{t,x) is small ? 

• What about an restricted three body problem in three dimensions? One 
obvious observation is that plane {x, y) is invariant for full 3D problem, 
hence we have symbolic dynamics also in a spatial problem. We would like 
to pose a more general question: Does there exists a symbolic dynamics for 
3D problem such the corresponding orbits are not all contained in the Sun- 
Jupiter plane ? Some preliminary numerical explorations in this direction 
can be found in paper |^ , 

• What about full 3-body problem? Will the symbolic dynamics established 
here continue to very small but nonzero mass of a comet? Some results in 
this direction for non-degenerate periodic orbits and generic bifurcations 
can be found in . 
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